Biophysics and Physicobiology

Streptavidin-biotin, avidin-biotin interactions are classical models for protein-ligand binding, yet the energetic changes accompanying biotin binding remain poorly resolved. Using fluorescent dyes as energy sensors, we show that biotin binding produces two distinct regimes occurring in parallel as concentration of biotin increases cooperativity and conformational rearrangements, wherein cooperativity is observed via exchange broadening of fluorescence linewidth and conformational rearrangements exclusively observed in emission energy. Where the first biotin binding creates the highest contribution to the emission energy. Further analysis of tetramer-tetramer only interactions revealed extremely long ranged intermolecular interactions extending to hundreds of nm. The intermolecular interactions become negligible only at concentrations of approximately 10 nM for both streptavidin and avidin. Affinity values estimated for these diluted samples were below 1 nM.

Amylose contributes to starch crystallinity, but the stability of packed amylose double helices in water at elevated temperature remains insufficiently characterized. Here, we used molecular dynamics simulations to test whether chain length affects the short-timescale stability of A-type amylose oligomers in water. Six systems differing in chain length (6, 12, or 24 glucose units per chain) and oligomer size (isolated double strand or dodecamer of six double strands) were simulated, and five independent 1-s production runs were analyzed for each simulated condition. Oligomers with six glucose units showed structural collapse accompanied by increased water penetration. By contrast, dodecamers with 12 or 24 glucose units largely retained packed double-helical organization over the simulated timescale, although fraying was observed at their ends. These results indicate that chain length and lateral packing strongly affect the early structural response of amylose-like crystalline segments in hot water. The present simulations do not establish the ultimate fate of longer oligomers at longer timescales, but they identify a relative stability difference that is relevant to molecular interpretations of hydration-driven disordering in starch.

Reconstructing time-resolved inter-residue distance distributions during protein functional dynamics in the solution state is known to be a difficult and important problem. This article presents a technique for extracting spin-spin (as a proxy for residue-residue) distance distributions on doubly-spin-labeled proteins from rapid-scan time-resolved Gd-Gd electron paramagnetic resonance (rs-TiGGER) spectra recorded near room temperature in solution at 240 GHz. We use a best-fit technique that convolves a dipolar kernel matrix with an intrinsic, non-dipolar-broadened (single-labeled) spectrum. The kernel incorporates the effect of solution-state tumbling on the dipolar broadening using a correlation function that bridges the static and rapidly tumbling regimes. We apply the technique to AsLOV2, a protein domain with a dark-state crystal structure that is well-known from X-ray crystallography, but a less well-characterized and disordered tertiary structure that manifests after photoactivation at 450 nm. Informed by principal component analysis, we assume that the underlying distance distribution may be approximated by a sum of two Gaussian distributions. The fits returned time-resolved, light-activated populations with mean distances of [Formula] (dark) and [Formula] (lit) in the wild type, and [Formula] (dark) and [Formula] (lit) in an N414Q mutant, with nearly complete unfolding (within fit uncertainty) of the active, light-sensitive fraction. The extracted distance distributions and their accompanying uncertainties are consistent within uncertainty with molecular dynamics simulations of the equilibrated protein structure.

Weak magnetic fields affect many biological processes across the tree of life, though the precise molecular sensors and pathways involved in such magnetoresponses remain mostly uncharacterized. Fluorescence is a useful tool for investigating magnetic field effects in flavoproteins, as their chromophores fluorescence intensity can be shown to depend on the spin states of electronic radical pairs. Here, we describe a four-state ordinary differential equation model to understand what parameter sets result in fluorescence contrast between spin states in photocycles with singlet and triplet radical pairs. We conclude that only certain sets of parameters result in the fluorescence intensity being a good proxy measurement for singlet yield. In particular, we observe that the illumination intensity required to obtain fluorescence contrast depends on the rate of the slow spin-independent radical termination reactions that recover ground-state oxidized fluorophores. Moreover, to observe a magnetic field effect in fluorescence intensity when an external magnetic field modulates the singlet yield, the illumination intensity must be strong enough such that photoexcitation is not the rate-limiting step. This understanding suggests that flavoproteins that do not exhibit magnetic field effects in their fluorescence emission under certain experimental setups may still be sensitive to weak magnetic fields in terms of function, as magnetosensitivity in fluorescence depends strongly on illumination conditions.

Heteromeric amino acid transporters (HATs) mediate essential amino acid flux across membranes, but the molecular dynamics of substrate translocation remain poorly defined for many family members. Here, using conventional and adaptive steered molecular dynamics (cMD and ASMD) simulations, we identify residue W230 in the b0,+AT transport channel as a dynamic gate that regulates arginine (Arg) influx through side chain flipping. By integrating dynamic network analysis with dynamical cross-correlation of residue motions, we show that regulatory signals propagate from the Arg binding site through transmembrane helix 5 (TM5), a connecting loop, and TM6 to reach W230. We propose a dynamic gating mechanism for b0,+AT - mediated amino acid transport. Arg binding at V186 triggers signal propagation that enhances cooperative interactions between W230 and Arg, driving the side chain flipping of W230. Our findings reveal a dynamic gating mechanism underlying b0,+AT - dependent Arg transport and suggest that residue-triggered side chain reorientation may represent a conserved and efficient strategy in transporter function. Author SummaryAmino acids are the essential building blocks of life, and their transport across cell membranes is vital for nutrition and cellular signaling. Heteromeric amino acid transporters (HATs) mediate this process, yet how they physically move substrates through the protein at the atomic level remains poorly understood. In this study, we used advanced computer simulations to observe, in unprecedented detail, how b0,+AT--a key HAT member--transports the amino acid arginine. Our simulations revealed that a single residue, tryptophan 230 (W230), functions as a molecular gate: its side chain flips open to allow arginine to pass and then closes behind it, ensuring one-way traffic into the cell. We further discovered that the initial binding of arginine sends a signal through specific structural elements (helices and loops) to trigger this gate opening. This work not only uncovers a dynamic gating mechanism for b0,+AT but also suggests that similar side-chain flipping events may represent a common and efficient strategy used by other transporters to control substrate movement. Our findings provide a new framework for understanding transporter function and could inform future drug design targeting these critical membrane proteins.

G protein-coupled receptors (GPCRs) are important targets for drug discovery because they are the largest and most diverse family of membrane proteins in the human body. As exemplified by {beta}2-adrenergic receptor ({beta}2AR), which is a typical class-A GPCRs, they are prone to denaturation and inactivation after solubilization. Although diverse techniques have been developed, current structural and functional analyses are still limited to proteins with relatively stable and high expression. However, some mutant variants and misfolded proteins involved in diseases have extremely low expression levels and may be difficult to analyze. To overcome these limitations, we established a novel analytical platform for evaluating the ligand-binding ability of full-length {beta}2AR as a model at the single-molecule level without purification and addressing challenges such as membrane proteins with low expression levels and structural instability. This method enables the direct use of unpurified receptors immediately after solubilization and allows us to use only a small amount of sample ([~]10 ng) for observation and to distinguish between specific and nonspecific ligand binding by fitting. Furthermore, we unveiled the physical properties of detergents on the structural stability of solubilized receptors and found that the lateral pressure within detergent micelles affects ligand-binding ability. Detergents that provided a fluid microenvironment were able to maintain ligand-binding ability for several days even after solubilization; conversely, detergents that provided a rigid microenvironment caused the protein to lose its activity earlier. Our method could be a promising tool for the structural and functional analysis of membrane proteins untargeted until now. Graphical abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=112 SRC="FIGDIR/small/740883v1_ufig1.gif" ALT="Figure 1"> View larger version (46K): org.highwire.dtl.DTLVardef@1ca0ed1org.highwire.dtl.DTLVardef@92e757org.highwire.dtl.DTLVardef@9223edorg.highwire.dtl.DTLVardef@74a100_HPS_FORMAT_FIGEXP M_FIG C_FIG

Interleukin-13 (IL-13) is an immunomodulatory cell signaling cytokine that has been implicated in neurodegenerative disease and chronic inflammation. IL-13 binds to its low and high affinity receptors, IL-13 receptor 1 (IL-13R1) and IL-13 receptor 2 (IL-13R2), respectively, with residence times that vary accordingly. As the binding kinetics of the cytokine-receptor complexes influence cellular responses, we employed the molecular dynamics (MD) simulation-based{tau} -random acceleration molecular dynamics method ({tau}RAMD) to compute relative residence times for wild-type (WT) IL-13 and 19 IL-13 mutants to the two receptors. Comparison with experimental kinetic data shows that the{tau} RAMD computations capture the trends in residence times. Analysis of simulated dissociation trajectories of the cytokine-receptor complexes reveals two distinct dissociation pathways of IL-13 from each of the receptors. This study thus pinpoints key determinants of the interaction of IL-13 with its receptors which could be targeted for therapeutic design. Statement of SignificanceCytokines are regulatory proteins that bind to cell surface receptors and thereby send signals to the cellular interior. Interleukin-13 (IL-13) is a cytokine that has a low and a high affinity receptor. It has important physiological roles, and its deregulation is involved in diseases such as atopic dermatitis and asthma. We employed a molecular dynamics simulation-based method to compute the effects of changes in the sequence of IL-13 on the lifetimes of complexes of IL-13 and its receptors. Comparison with experiments supports the validity of the computational approach and analysis of the simulations reveals two distinct ways in which IL-13 dissociates from each receptor. These results thus provide a map for targeting IL-13 - receptor interactions for the design of therapeutics.

The ability of small molecules to cross the blood-brain barrier (BBB) remains a major bottleneck in neurotherapeutic development. While experimental assays and machine learning approaches provide approximate permeability estimates, they lack atomistic insight into the underlying transport mechanisms. Here, we employ all-atom molecular dynamics simulations of a compositionally realistic BBB lipid bilayer to characterize the passive permeation of two bioactive propolis-derived compounds, Caffeic Acid Phenethyl Ester (CAPE) and Artepillin-C (ARC). Using steered molecular dynamics and umbrella sampling, we computed free energy profiles, diffusion coefficients, and permeability metrics across the membrane. CAPE encounters a modest barrier at the lipid headgroup region but minimal resistance within the hydrophobic core, resulting in a low free energy barrier ([~]2-3 kcal/mol) and favorable permeability (logP_eff {approx} 0.28). In contrast, ARC exhibits a substantial energetic barrier within the membrane core, leading to high resistivity and strongly unfavorable permeability (logP_eff {approx} -10.91). The heterogeneous lipid model reproduces experimentally consistent membrane properties and reveals how lipid composition modulates transport energetics. These findings provide mechanistic insight into BBB permeability and demonstrate the utility of atomistic simulations for guiding the design of neuroactive therapeutics.

Magnetic treatment in tomato seeds (Solanum lycopersicum L.) has been studied as a biotechnological technique to induce a reduction in germination times and enhance plant development. However, the modified cellular mechanisms involved in the reduction of germination times or the improvement of development parameters are not yet clearly established. To explore a possible altered cellular mechanism, the effect of homogeneous static magnetic fields on the structure of the cyclic nucleotide-gated channel 6 (CNGC6), the modification in the organization of POPC lipids in the plasma membrane, and changes in calcium ion mobility were evaluated. For this purpose, coarse-grained molecular dynamics simulations were performed using the Martini 3 model in GROMACS, applying five different magnetic flux densities (0.000, 0.001, 0.010, 0.100, 1.000, and 10.000) T over 1 000 ns. The results showed an anisotropic effect in the longitudinal direction of the protein, which generated heterogeneous behavior among the chains of the homotetramer; this altered the conformation of the CNGC6 channel and modified the pore bottleneck. In contrast, no significant changes were observed in the conformational order of the POPC phospholipid chains. As a preliminary, single-replicate exploratory study, these results suggest that homogeneous static magnetic fields may induce specific structural modifications in the CNGC6 ion channel of Solanum lycopersicum L. without compromising the integrity of the lipid bilayer or the dynamics of ion transport within the analyzed timescale; these preliminary findings provide a molecular-level structural basis for future experimental and computational investigations of magnetic field effects on plant cyclic nucleotide-gated channels.

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Attenuated total reflectance Fourier-transform infrared (ATR-FTIR) spectroscopy of proteins in aqueous solution is often limited by water absorption and other optical artifacts. To overcome these limitations, we evaluated the structural features and hydrogen-deuterium exchange (HDX) kinetics of the -helical protein GCN4 in both hydrated (wet) and vacuum-dried (dry) states. While solvent heavily mask the second-derivative spectra of wet samples, vacuum drying yielded a thin, protein-rich film on the ATR crystal, significantly enhancing the signal-to-noise ratio and resolving the protein features without altering the native structure. Dry-state analysis clearly resolved the Amide I, Amide II, and deuterium-shifted Amide II' (1450 cm-1) bands. Notably, second-derivative analysis of the dry spectra of the HDX samples revealed a bimodal Amide I distribution consisting of a stationary band at 1653 cm-1 from the solvent-inaccessible regions and an isotopically sensitive band shifting from 1648 cm-1 to 1644 cm-1 from solvent-accessible regions. These results demonstrate that vacuum-dried ATR-FTIR spectroscopy effectively eliminates solvent masking, providing the spectral clarity required to resolve discrete -helical sub-populations after deuteration.

The Gene Version Iteration Hypothesis (GVIH) proposes that mutant genes may originate from the Y chromosome, traverse through the X chromosome to autosomes, undergo interchromosomal transfer, and potentially return to the Y chromosome via the X chromosome. This hypothetical closed transmission loop may facilitate the storage, screening, and elimination of different versions of mutant genes. The hypothesis comprises five core propositions: (1) Mutation reservoir: The Y chromosome may serve as a specialized carrier for generating mutant genes, characterized by elevated mutation rates, reduced gene density, and accelerated evolutionary dynamics; (2) Closed-loop transmission: Mutant genes may follow a unidirectional pathway Y[->]X[->]autosomes[->]X[->]Y, forming a complete transmission circuit; (3) Coexistence of multiple versions: A single functional gene may exist in multiple versions across different chromosomes, constituting a dynamic gene version library; (4) Reproductive screening: Environmentally adaptive gene versions may persist across generations and potentially migrate to upstream chromosomes, while maladaptive versions may be eliminated; (5) Terminal elimination: Gene versions reaching the Y chromosome may undergo elimination processes, potentially preventing version monopolization and maintaining evolutionary dynamics. This hypothesis provides a novel framework for understanding adaptive evolution at the genetic level. If empirically validated, it may offer new insights into the molecular mechanisms underlying certain genetic phenomena and evolutionary processes.

Members of the transforming growth factor-{beta} (TGF-{beta}) superfamily, including inhibins and activins, are structurally related glycoprotein dimers that regulate reproductive and endocrine signaling. Their high degree of molecular similarity presents challenges for label-free analytical discrimination. To evaluate the ability of Raman spectroscopy to distinguish closely related TGF-{beta} superfamily proteins based on intrinsic vibrational fingerprints. Raman spectra of recombinant human Inhibin -subunit, Inhibin B ({beta}B homodimer), and Activin A ({beta}A--{beta}A) were acquired using confocal Raman microscopy with 532 nm excitation. Spectra were baseline-corrected, area-normalized, and analysed using principal component analysis (PCA). Distinct spectral signatures were observed across the 500--1800 cm-1 region. Differences within the S--S stretching region (500--550 cm-1) were consistent with variations in disulfide-bond environments, with the Inhibin -subunit exhibiting the highest relative intensity in this region. Variations in the amide I band (1600--1700 cm-1) suggested differences in protein secondary structure, while aromatic amino acid vibrations provided additional discriminatory features. PCA revealed clear clustering and separation of all three protein classes based on their Raman fingerprints. Raman spectroscopy enables label-free differentiation of structurally related endocrine glycoproteins and demonstrates potential for the structural characterization and classification of inhibin and activin proteins within the TGF-{beta} superfamily.

Glucocorticoids are among the most widely prescribed drugs globally due to their potent anti-inflammatory and immunosuppressive actions. These effects are primarily mediated by the glucocorticoid receptor (GR), a ligand-activated transcription factor that translocates from the cytoplasm to the nucleus to regulate hundreds of genes. Although nuclear entry is a prerequisite for its genomic response, the mechanisms governing this process remain unresolved; specifically, whether the receptor translocates as a monomer or a dimer remains a subject of significant controversy. Here, we employed the pair correlation function (pCF) approach to quantify the nuclear translocation of single fluorescent GR molecules in live cells. This minimally invasive method identifies correlations between intensity fluctuations generated by molecules moving from the cytoplasm into the nucleus. Our results demonstrate that GRs quaternary structure and conformation modulate GR transport. While GR monomers rely exclusively on passive diffusion, GR dimers also utilize the microtubule-dynein machinery for active transport, proving that dimerization can precede nuclear import. Furthermore, the perinuclear vimentin network facilitates faster translocation by constraining actively transported dimers near nuclear pores. Collectively, our work reconciles contradicting reports regarding GR stoichiometry during import by demonstrating that both monomers and dimers translocate, albeit through distinct mechanisms. Importantly, these results reopen the door for a microtubule-dependent, heterocomplex-independent model of GR translocation, suggesting that the cytoskeleton is an integral, yet overlooked, component of the GR signaling pathway.

We demonstrate the use of fluorescence correlation spectroscopy (FCS) to characterize fluorescently-labeled protein production. We use cell-free protein synthesis to express the protein YFP-CopB, a fusion of Chlamydia Outer Protein (Cop) B and Yellow Fluorescent Protein (YFP). CopB is a [~]50 kDa protein believed to have a critical role in chlamydial infection.1 After adding a plasmid encoding YFP-CopB to an E. coli cell-free lysate, protein expression begins. We track the cell-free reaction over several hours using the EI-FLEX, a commercial instrument with FCS capability. As protein is expressed over time, YFP-CopB increases in concentration, and the EI-FLEX detects an increase in fluorescent signal above the background of the cell-free lysate. The FCS data collected gives information about the size, aggregation tendencies, rates of production and fluorescent protein maturation, and concentration of the YFP-CopB produced. The use of FCS concurrent with cell-free synthesis presents a simple method to characterize proteins of interest as they are produced without the need for purification.

The process of drug unbinding is of immense importance in the field of biophysics and therapeutics. The behavior of these systems is greatly influenced by their thermodynamic and kinetic properties. Therefore, it is crucial to accurately estimate the ligand binding free energies and rate of ligand dissociation, yet these processes are often governed by rare event transitions that lie beyond the reach of standard brute-force molecular dynamics simulations. While enhanced sampling simulations offer a solution, their efficacy is strictly contingent upon the selection of appropriate collective variables (CVs) which is non-trivial for complex systems like protein-ligand complexes. In this study, we present a method to derive optimized CV from transition state region (TS) via an interpretable machine learning (ML) model, Elastic Net. By employing some physically intuitive order parameters, the derived optimized CV from the TS-region greatly accelerate ligand binding-unbinding transitions and achieves rapid free energy surface (FES) convergence across diverse systems including buried and solvent exposed active sites such as Trpsin-benzamidine complex, host-guest systems and sodium epoxidase etc. Intriguingly significant contribution of the ligand hydration is found in the optimized CV which depicts crucial role of solvent in driving ligand binding-unbinding transitions. The estimated binding free energies for different protein-ligand complexes match quite well with experiments, while maintaining a low computational cost. The derived optimized CV is also used to calculate the ligand residence times across different systems and calculated residence times are within the experimental range for all systems, again with very little computational costs. Moreover, we show that the optimized CV constructed from TS region via an interpretable ML model is transferable across diverse systems, offering a robust and scalable framework for drug discovery and investigation of complex biomolecular recognition.

Far-Red (FR) Light Photoacclimation (FaRLiP) enables cyanobacteria to extend photosynthetic activity into the far-red region by extensively remodelling Photosystem I (PSI), including the replacement of several core subunits with paralogs that coordinate the red-shifted chlorophyll f (Chl f). The binding positions of Chls f are still a matter of debate, with the most recent structural findings supporting the location of a single Chl f molecule within the reaction centre (RC) at the so-called A-1B site. This was in turn suggested to strongly affect electron transfer (ET) directionality leading to an almost monodirectional transfer along the B branch in FR-PSI RC. Here, we directly probe ET in FR-PSI by characterising the photogenerated [P700A1-] spin-correlated radical pair using complementary pulse and Time-Resolved (TR) Electron Paramagnetic Resonance (EPR) spectroscopy at cryogenic temperature. Electron spin-echo decay kinetics are distinctly biexponential, indicating the formation of two charge-separated states. Consistently, out-of-phase ESEEM traces are quantitatively described by two modulation frequencies arising from different dipolar interactions, while TR-EPR spectra are accurately simulated by the combined contributions of [P700A1A-] and [P700A1B-] radical pairs. These results provide direct spectroscopic evidence that both the A and B branches remain photochemically active in FR-PSI. The conservation of bidirectional ET, even when considering the presence of a single Chl f molecule in the RC, further implies that the two radical pairs originate from a common primary electron donor. This finding identifies P700 as the most likely primary donor and argues against a mechanism in which the RC Chl f initiates charge separation.

A universal mechanism in molecular evolution is functional and structural divergence of members of a protein family. The ability of AlphaFold to predict atomic-resolution protein structures promises to accelerate insights into this process. We study the interplay of changes in sequence, structure, and function in photoactive yellow protein (PYP), a family of bacterial blue light photoreceptors. Halorhodospira halophila contains two PYP homologs that diverged to 60% sequence identity, differ 100-fold in the lifetime ({tau}pB) of their pB signaling intermediate, and display altered peak wavelengths ({lambda}max) for color sensing. We resurrected ancestral PYPs and determined these properties along the resulting recapitulating evolutionary divergence. The resurrected ancestral PYP is functionally similar to PYP1, indicating divergence on the path to PYP2. AlphaFold predictions for PYP2 and these ancestral proteins revealed the absence of structural changes compared to the crystal structure of PYP1. To experimentally validate these predictions, we optimized second-derivative Fourier transform infrared (FTIR) spectroscopy. The FTIR spectra of PYP1 and 2 and their resurrected ancestral proteins demonstrated clear differences in their secondary structure. These results demonstrate an important limitation of AlphaFold and show how ancestral sequence reconstruction combined with spectroscopic approaches yields insights into divergence in a protein family.

Phosphatase and tensin homolog (PTEN)-induced putative kinase 1 (PINK1), a key regulator of mitophagy, has been linked to the pathogenesis of Parkinson's disease (PD). PINK1 recruits Parkin, an E3 ubiquitin ligase, triggering mitophagy in response to mitochondrial damage. During mitophagy, the quantity, stability, and activity of PINK1 must be strictly regulated; however, the mechanisms governing these parameters under cellular stress are still unclear. Herein, we determined the structural basis for PINK1 maturation mediated by heat shock protein 90/cell division cycle 37/FK506-binding protein 51 (HSP90/CDC37/FKBP51) chaperone complex. We identified PINK1-associated proteins using liquid chromatography-tandem mass spectrometry (LC-MS/MS) and determined the structures of the complexes using Cryo-Electron Microscopy (Cryo-EM). Results showed that FKBP51 potentially interacts with a conserved leucine-proline-phenylalanine (LPF) motif on the activation loop of PINK1 and negatively regulates PINK1 functions in mitophagy. A PINK1 mutation located at the FKBP51 recognition site is linked to mitophagy deficiency, which can be partially rescued by specific inhibition of FKBP51. These findings reveal a general mechanism for PINK1 recognition by the HSP90/CDC37/FKBP51 chaperone complex and suggest a potential approach for upregulating PINK1 activity, which is impaired in PD.

N-acetylaspartate (NAA) is the most abundant neuron-enriched acetylated metabolite in the mammalian brain, but its metabolic purpose remains unresolved. We developed a simplified kinetic model of mitochondrial aspartate metabolism to test whether NAA synthesis by aspartate N-acetyltransferase (ASPNAT) acts as a thermodynamic relief valve for mitochondrial aspartate aminotransferase (AAT) under the low-oxaloacetate (OAA) conditions expected in neuronal mitochondria. In the mitochondrial-compartment model, ASPNAT lowered steady-state mitochondrial aspartate from 141 to 105 M and increased net forward AAT flux by 30.9%. The relative AAT-relief effect was largest when OAA and aspartate-glutamate carrier 1 (AGC1/Aralar1)-mediated export were both low, whereas acetyl-CoA availability controlled the substrate-supported capacity for NAA synthesis. That places the relief effect in a narrow regime where product removal matters most. ASPNAT titration produced a graded, concentration-dependent response rather than a binary on/off response. Energetic comparisons showed that the gain in AAT-linked support comes at a modest acetyl-CoA cost, which makes NAA synthesis easier to sustain in carbon-replete states than in carbon-poor ones. Some studies have suggested a secondary cytoplasmic site of NAA synthesis, and we therefore examined how the network response changed with a change in ASPNAT topology. Mitochondrial matrix ASPNAT increased forward AAT flux by 53.32%, whereas cytoplasmic ASPNAT decreased ASPNAT flux by 17.8%. Allowing OAA to vary preserved the positive ASPNAT-dependent relief of AAT flux, but because this simplified extension produced unrealistically low absolute fluxes, it is interpreted as a robustness check on the direction of the mechanism rather than as a prediction of physiological metabolic rates. These results identify mitochondrial NAA synthesis as a plausible thermodynamic relief valve for mitochondrial AAT and define a directional prediction that could test whether severe metabolic stress reroutes effective ASPNAT-linked aspartate metabolism.